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Abstract. We present an analysis of the (5-Scuti star object HD 174936 (ID 7613) observed by CoRoT during the first short 
run SRcOl (27 days). A total number of 422 frequencies we are extracted from the light curve using standard prewhitening 
techniques. This number of frequencies was obtained by considering a spectral significance limit of sig=10 using the software 
package SigSpec. Our analysis of the oscillation frequency spectrum reveals a spacing periodicity of around 52 yuHz. Although 
modes considered here are not in the asymptotic regime, a comparison with stellar models confirms that this signature may 
stem from a quasi-periodic pattern similar to the so-called large separation in solar-like stars. 
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1. Introduction 

The (J-Scuti stars are intermediate-mass pulsating variables 
with spectral types ranging from A2 to FO. They are located on 
and just off the main sequence in the lower part of the Cepheid 
instability strip (luminosity classes IV & V). The (5-Scuti stars 
are considered as particularly suitable for asteroseismic studies 
of poorly known hydrodynamical processes occurring in stellar 
interiors, such as convective overshoot, mixing of chemical el- 
ements, and redistribution of angular momentum (Zahn 1992), 
to mention only the most important ones. Because of the com- 
plexity of the oscillation spectra, their pulsating behaviour has 
not bee n understoo d very well. The reader is referred to the re- 
view bv iCoxl(l2002h for a more detailed description of unsolved 
problems in stellar pulsation physics. 

The 5-Scuti star HD 174936 was uninterruptedly observed 
by CoRoT during 27 days of observation. The very precise 
space photometry supplied by the CoRoT mission gives us the 
possibility of working with a lower limit in amplitudes of ~ 14 
//mag. As a consequence of such precision, never reached by 
any ground-based observation campaign, we obtained a wide 
range of oscillation frequencies ([0.05, 100] cd"'). This and 
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the huge number of detected frequencies for such stars will help 
us in understanding the internal structure of intermediate-mass, 
main-sequence stars. 

One of the aims of the mission is that such high resolution 
can e nable us to detect some periodicities like rotational split- 
ting dBaglin et al.ll2006bl) . Using the CoRoT achievements, we 
can carry out statistical studies and look for other possible peri- 
odicities that could help us to better understand the oscillation 
spectra observed. 

The paper is structured as follows: an explanation of the 
data analysis and the frequency content is described in Sect. 
121 a description of the method used for identifying the peri- 
odic spacing distribution of the extracted frequencies is given 
in Sect. [3j equilibrium models and oscillation computations to 
check the origin of the signature noticed in the observations are 
described in Sect.|4l a brief discussion is carried out in Sect.|5l 
and some conclusions are reported in Sect.|6l 

2. Data analysis 

The field (5-Scuti star HD 174936 was observe d during the first 
30-day short run SRcOl jBaglin et al.ll2006iJ) of CoRoT. The 
final time span of the collected dataset was AT=27.2 days, 
with a sampling of one point each 32 sec. This means about 
73440 datapoints as the expected timeseries. However, the fi- 
nal dataset consisted of 66057 datapoints after removing those 
points considered unreliable. With this, the Rayleigh frequency 
resolution is (1/A T)=0.037 cd"', and an oversampling of 20 
should correspond to a frequency spacing of 0.0018 cd This 
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Fig. 1. Spectral window of the short-run CoRoT dataset on 
HD 174936. 
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Fig. 2. The 422 frequencies exctracted for HD 174936. 



is equivalent to the frequen cy spacing in mode " High" in the 
program package Period04 ( Lenz & Bregerll2005 ). 

First, the data were detrended dAuvergne et al. IL2009') per- 
forming a linear fit to the light curves. Second, Period04 was 
used for a preliminary inspection of the periodograms. Later, 
the timeseries was analysed using the computer program pack- 
age SigSpec (Ree gen 2007) . 

Nevertheless, the timeseries was previously investigated 
with Period04 for the first 20 peaks. They were identical to 
those achieved with SigSpec. This test has also been carried 
out successfully for the CoRoT target ID 123 (HD 50844, 
Poretti et al.l l2009') with AT=56.7 days and 140016 datapoints. 
In that case, Period04 was used to investigate the first 200 
peaks. The agreement between the two methods was successful 
in 99% of the cases. 

Figure [T] shows the spectral window of HD 174936, which 
is typical of all the targets observed by the CoRoT satellite in 
the same run. As can be seen, the figure does not show the typ- 
ical aliases at 1 cd"' or the power levels that are common for 
ground-based data. On the contrary, all the aliases are related 
to effects produced by the satellite and its orbital frequency 
(fj=13.972 cd"'), and their power levels are much lower than 
those usual for ground-based data. The alias peak at about 2 
cd"' comes from the South Atlantic Anomaly (SAA) crossing, 
which occurred twice each sidereal day. This means the alias 
peaks at 2.005 and 4.01 1 cd"', respectively. 

To avoid problems with power close to zero frequency, the 
analysis with SigSpec was carried out in the range 0.05-100 
cd~^ In the case of Period04, the limit commonly used to con- 
sider one datapoint as significant is the amplitude signal-to- 
noise ratio S/N>4.0. In the case of SigSpec, the parameter used 
for significance is "sig" (- spectral significance), and the limit 
used as default is sig=5.0. This should be equivalent to about 
S/N=3.8 (an d sig=5.46 should be approximately equivalent 



to S/N=4.0) jReegenl 120071; iKallinger et alj l2008h . However 



we used a much more conservative limit, and our calculations 
with SigSpec stopped when sig=10.0 because the correspond- 
ing S/N values, determined using Period04 on the residuals, 
were much lower than expected. This is probably caused by a 
high number of peaks still remaining among the residuals in 
the region of interest. This was explained in much more detail 



in similar recent works for other CoRoT targets (HD 50844, 
Poretti et al ]l2009l HD 49434, Rodriguez et al., in preparation). 

The limit of sig= 10 was achieved after removing 422 peaks. 
This means a level of about 14 ppm for the smallest amplitudes. 
The first eleven frequencies, are listed in Table [T] along with 
the most relevant parameters. Column 6 lists the S/N values 
corresponding to each peak. These S/N values were calculated 
using Period04 on the residual file provided by SigSpec. Each 
S/N value was calculated on a box of width=5 cd"' centred 
on the corresponding pe ak, as is usual for this type of variable 
(IRodriguez et al ] l2006allbL . 

The range of statistically significant detected frequencies 
goes from some value close to zero up to about 70 cd ', i.e. 
800/iHz (1 cd ' = 11.57;uHz). But those with the highest am- 
plitudes are grouped around 35 cd ' (400/iHz, see Fig.|2]l. In 
the following theorical analysis, we work in /iHz. 



3. Looking for periodic patterns in the extracted 
frequency set 

Mode identification, mainly for (5-Scuti stars, is the present bot- 
tleneck for the progress in asteroseismology. Existing mode 
identification techniques are based on spectroscopy or multi- 
colour photometry not available in CoRoT data. Considering 
that we have unprecedented, rich spectra, we seek to reveal 
characteristic structures in the frequencies distribution and to 
start with quasi-regular spacings. 

We consider the frequency set as a series of Dirac's (5's of 
equal amplitude. Then, the technique is simply to use the list 
of frequencies with unit amplitude and to normalise its Fourier 
transform to unity at zero frequency spacing. If we have a given 
periodicity in the frequency set {6f), then its Fourier transform 
will show a clear periodic structure on the corresponding in- 
verse scale. 

We started this analysis assuming the hypothesis that, se- 
lecting a subset with the highest peaks of the entire frequency 
set, we are actually considering mainly the modes of the lowest 
{ values. The visibility of the modes decreases approximately 
as {^^-^ or depe nding whether the { degree is odd or even 
( Dziembowskill 1977 ). We analysed several subsets with an ar- 
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tification 


Frequency (yjHz) 


Amplitude (mmag) 


Phase (rod) 


Spectral significance 


S/N 


Fl 


377.30 


2.12 


2.40 


7225.844 


659.008 


F2 


412.71 


1.02 


2.49 


3624.781 


336.581 


F3 


414.62 


0.72 


-0.41 


2155.327 


236.220 


F4 


360.08 


0.57 


2.67 


1784.588 


174.791 


F5 


339.22 


0.55 


1.18 


1714.847 


169.363 


F6 


367.96 


0.53 


-0.82 


1862.847 


164.484 


F7 


321.92 


0.36 


-1.36 


908.218 


111.772 


F8 


387.63 


0.32 


-2.55 


825.756 


103.085 


F9 


359.47 


0.30 


2.23 


713.637 


92.766 


FIO 


385.89 


0.29 


-0.05 


725.392 


93.201 


Fll 


158.01 


0.28 


-2.44 


707.305 


85.501 
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Fig. 3. Power spectrum for various subsets of frequencies, selected by amplitude. The solid line represents the power spectrum 
when the highest 50 frequencies are selected; the dashed line corresponds to the highest 100 frequencies; the dotted line, for 200; 
and dot-dashed hne for aU of them. The peak corresponding to the large separation (25.8 //Hz) and its sub-multiples are labelled 
(see text). 



bitrary number of frequencies but always including the highest 
peaks. In Fig.[3]the results are given for different sets of peaks. 
It can be seen that the more frequencies in the subset, the flatter 
the peaks, although the local maxima remain at the same posi- 
tion. This indicate that, when we only select the highest peaks. 



we are enhancing some periodicities corresponding to a few t 
values. While, adding frequencies in the analysis, periodicities 
of all modes are powered and no pattern can be significantly 
distinguished from others. 
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When we select 50 frequencies, a periodic pattern can be 
recognised. Furthermore, the clear signature of a Dirac's comb 
is also present at 25.8, 13, 8, 6.5, and 4.9 juHz. That indicates 
that we probably have a Dirac's comb of 25.8 yuHz, which is the 
lowest divisible number The appearance of these peaks can be 
a clue to the confirmation that they are being caused by a pe- 
riodic structure in the frequency set. Other peaks around 25.8 
fiHz can be explained by the separation not being strictly pe- 
riodic but quasi-periodic, as can be expected for these modes 
(see Sect. nil. This can be the imprint of a large separation-like 
structure around 52 /iHzQ 

On the other hand, the rotational velocity measured for this 
star is high (vsin/ ~ 170 km/s, see Sect. |4]i giving a mini- 
mum rotational splitting of ~ 20 yuHz. It is then possible that 
the observed frequency pattern does not come from the large 
separation but from the rotational splitting. The pattern due to 
rotational effects should be observed only with { > I and a sta- 
tistically significant number of ;« i^Q values. However, because 
we mainly select the lowest { values, the large separation is 
more probable than the rotational splitting. In the next section, 
we discuss the validity of this hypothesis using a theoretical 
model representative of the star 

The other peaks in Fig. [3] probably come from periodic 
structures produced by the rotation and/or the small separation, 
mainly occupying the lowest part of the power spectrum. Works 
are in progress for understanding the meaning of these peaks. 



4. Periodic structures in theoretical spectra 

The periodic pattern structures are familiar in solar-like oscil- 
lators, which have reached the asymptotic regime, but not in 6- 
Scuti stars, which pulsate around the fundamental radial mode. 
To understand the physical origin of the present observed pe- 
riodic pattern shown in the previous section, we computed a 
representative model of this star. 



4.1. Physical parameters 

We consider the physical parameters of the star HD 174936, 
which are listed in Ta ble [2] These va l ues ar e taken from 

Teff, logg, 

photometry 



Table 2. Physical parameters of t he star with its uncert ainties 
to construct the theoretical model dCharpinet et al.ll2006 ). 



"CorotSky Database" dCharpinet et all l2006h . 
and \Fe/H] were derived from Stromgren 



( Hauck & Mermilliod 1998h . while the rotational velocity 
(vsin; = 169.7km/ s) was determined from high-resolution 
spectroscopy. The spectrum was taken in June 2, 2004 using 
the FEROS spectrograph attached to the ESO 2.2 m-telescope 
at La Silla (Chile). It was obtained in the framework of the 
mission prep aration and is available at the GAUDI archive 
( Solano et al.ll2005.) . On the other hand, it has been shown that 
rapid rotation should be considered when calculating the stel- 
lar physical parameters from photometry. IMichel et al.l (119991) 
propose a method for determining the effects of rotation and 
geometry (angle of inclination) on photometric parameters for 
stars in clusters. For i5-Scuti stars, that method was further de- 



Star ID 


Teff (K) 




[Fe/H] 


vsin; (km/s) 


HD 174936 


8000±200 


4.08±0.2 


-0.32+0.2 


169.7 



veloped by lPerez Hernandez et al. I (ll999h . showing that uncer- 
tainties of around 100-150 K in effective temperature, and ~ 
0.10 dex in logg, can be found fo r moderately rotatin g stars. 
That result was later confirmed bv lSuarez et al.l (120021) . In the 
present case, considering the absence of additional information 
on the inclination angle of the star, uncertainty boxes of 200 K 
in the HR diagram in Te/f and ~ 0.2 dex in logg were adopted. 
An uncertainty of ~ 0.2 in metallicity was taken. 

This is a main-sequence intermediate-mass star (~1.6 Mq). 
The hydrogen abundance in t he core is Xg ~ 0.3. But it is 



younger than the star studied in lPoretti et al.l (l2009h . 



4.2. Theoretical models 

Using the physical parameters measured for the star and given 
in the previous subsection as input, we computed a theoreti- 
cal model represe ntative of the star. The evo lutionary code ce- 
SAM (.MoreL1997l:lMorel &Lebretonll2008l) and the pu lsation 



' Frequency values for ^ = 1 modes are placed around the centre of 
two consecutive f = modes, so we expect to find a minimun period- 
icity of the half value of the large separation. 



code GRAco dMoya et all 120041: iMoya & Garridol l2008h were 
used as numerical codes to calculate frequencies, growth rates 
and other physical quantities. 

For the model, we computed the oscillation spectrum (ra- 
dial and non-radial), from { = to 3, in the range of observed 
frequencies. Then, we calculated the large separation defined 
as Av - v„+i/ - v„j. In Fig.|4]the large separation calculated 
this way is depicted. It is quite evident that the large separation 
does not vary much but is bounded within a frequency range 
of 10 juHz. In this particular case, the two discrepant points be- 
longing to f = 2 correspond to an avoided crossing. Similar 
results are expected for other higher i values. 

However, this is not critical for the detection of regulari- 
ties using Fourier transform due to its statistical character. In 
Fig.|5]we show an identical analysis to the one for the observed 
frequencies (see previous section). It can be seen that we have 
similar results for the theoretical frequencies for a given repre- 
sentative model (even in a nonrotating model). This indicates 
that such regularities, even in the non-asymptotic regime, may 
be discerned by the Fourier transform technique explained in 
Sect. [3] This opens a new window in the analysis of 5-Scuti 
stars. 

5. Discussion 

We have shown that a periodic pattern, possibly correspond- 
ing to large separation, is found in the frequency spectrum 
observed for the star (see Sect. O. We have coiToborated it 
through a theoretical model (see previous section). Clearly, 
the large separation found in the studied regime does not have 
the same physical meaning as found in the asymptotic regime 
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Fig. 4. Frequency dependance of the large separation for the 
three lowest values of (: 0, 1, and 2. Solid line is for ( =0, 
dashed line for t - I, and dot-dashed line for { -2. Large sepa- 
ration is then estimated using a Fourier transform method (see 
Sect.EJ. 



Fig. 6. HR diagram containing two uncertainty boxes corre- 
sponding to the Icr extremes of the metallicity ([Fe/H] = - 
0.12 and -0.52), with two evolutionary tracks. They correspond 
to the central model (M = 1.63 M©) at the central metallicity 
([Fe/H] = -0.32), one with rotation and the other without. 



Power spectrum of the frequency ([100,700]^iHz) distribution 




Periodicities in tlie frequency distribution (tlHz) 



Fig. 5. Power spectrum of the oscillation frequencies of a 
model obtained using the technique described in Sect. |3] We 
use only the range of observed data set, i.e., up to 700 yuHz. A 
clear peak (and sub-multiples) is observed at 26.61 //Hz, which 
corresponds to a large separation of ~53 jt/Hz. 



dChristensen- Dalsgaardl2003l) . ' iowever, it is dependent on the 
internal structure of the star, and we can study how this sig- 
nature is robust and sensitive to evolution or mass, in order to 
open up new perspectives about its diagnostic power. 

To do so, we constructed seismic models that lie within the 
uncertainties of the physical parameters measured for this star 
That is, we computed models to reproduce the HR locations 
corresponding to the four corners (1 it) and the centre of the 
uncertainty box. After that, for each equilibrium model, we 
computed its oscillations. This procedure is followed for mod- 
els built with the observed metallicity and the two maximum 



errors (1 cr) considered. Then, we reduced the corners of the 
uncertainty box to half and computed the corresponding mod- 
els. Finally, we calculated the large separation for all them. The 
uncertainty boxes (1 cr) for the 1 cr extremes of the metallicity 
are depicted in Fig. |6] Two evolutionary tracks of the central 
model (M = 1.63 M©) with and without rotation are plotted, 
too. 

We can see in Table [3] the results of all these calculations. 
The large separation is close to the value found in the ob- 
served frequency set for the central models at each metallicity. 
Typically, it is higher for less massive models, and when reduc- 
ing the half -bottom area of the uncertainty box it varies faster 
than the half-upper one. 

On the other hand, we computed the instability range for 
these models (see Fig. |7]). Some models covers the range 
(around 400 juHz) of the highest amplitudes in the frequency 
set, others do not. But not one of them gives the necessary in- 
stability range to represent the full range in the observations. 
It is assumed that observed frequencies ([0.5, 900] juHz), at 
least in ground-based data, must be within the instability range. 
However, as far as we know, no theory allows such a huge range 
of excited frequencies, even including the time-dependent con- 
vection theory. This is a new challenge for theory that should 
be improved in order to predict the instability range observed 
in the new era of space telescopes. 

As seen, the observed frequency set contains a lot of infor- 
mation (small separation, rotational splitting, etc.), and work 
is progressing. We hope that such a large number of frequen- 
cies will allow us to obtain precise models representative of 
this star and to improve our understanding of the interiors of 
intermediate-mass stars. 
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Table 3. Models of the uncertainty boxes for each metallicity calculated with mass and large separation. 





First box 






Second box 




[FelH] 


M(Mo) 


Ay (^Hz) 




M(Mo) 


Av (p.Hz) 




1.63 


53 




1.63 


53 




1.89 


37 




1.77 


44 


-0.32 


1.77 


38 


-0.32 


1.71 


44 




1.51 


82 




1.58 


60 




1.43 


71 




1.54 


58 




1 ^ 1 






1 s 1 
1 .J 1 






1.76 


39 




1.64 


44 


-0.52 


1.62 


41 


-0.52 


1.58 


45 




1.42 


74 




1.48 


60 




1.33 


73 




1.43 


60 




1.76 


51 




1.76 


51 




2.05 


37 




1.92 


43 


-0.12 


1.90 


39 


-0.12 


1.85 


43 




1.66 


71 




1.73 


58 




1.55 


73 




1.67 


58 
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Fig. 7. Frequency versus growth rate and amplitude for the 1 .63 
M0 model. Growth rate (solid line, left axis) shows the stability 
range (< 0), and one vertical line of the length of the amplitude 
for each frequency is represented (right axis). 



6. Summary & conclusions 

This work has presented an analysis of the 5-Scuti star ob- 
ject HD 174936 (ID 7613), observed by CoRoT during the first 
short run SRcOl (27 days). A total number of 422 frequencies 
(which reach up 700 juHz) are extracted from the light curve 
using standard prewhitening techniques. This represents one of 
the largest number of detected frequencies for a i5-Scuti ever 
obtained so far. 

We combined the classical seismic analysis with the use of 
statistical properties of the modes distribution in the observed 



frequency range. In particular, we have found periodic patterns 
in the observed frequency spectrum, which were not expected 
for this kind of pulsating star. We find a frequency distribution 
of the modes that seems to correspond to the large separation, 
which is about 52 juHz. This result is supported by an equiva- 
lent analysis performed with a representative model of the star 
All these results provide new prospects for the asteroseis- 
mology of (5-Scuti stars. In particular, the large separation 
found for this star opens the possibility of performing simi- 
lar studies to those for solar-like stars. Although we have not 
explored the low-frequency region of Fig. [3] we cannot dis- 
card the possibility that relevant information concerning rota- 
tion and small differences are hidden in this peak forest. 
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